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ABSTRACT 



The United States Air Force uses a nonlinear programming model to assess the 
utilization of weapons and sorties needed to achieve a maximum value of destroyed 
targets in a multi-period, Theater- Level conflict. The current model is modified by con- 
straining the consumption of weapons. Alternate objective functions are introduced. 
Their meaning and influence on the optimization is compared. An increase in the worth 
of destroyed targets is gained if the model can more flexibly utilize weapons than is 
currently the case. The optimization can be further improved if all time periods are 
considered simultaneously while assigning sorties to targets, rather than the current 
myopic approach. 



THESIS DISCLAIMER 



The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within the 
time available, to ensure that the programs are free of computational and logic errors 
they cannot be considered validated. Any application of these programs without addi- 
tional verification is at the risk of the user. 
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I. INTRODUCTION 



In 1988 the United States Air Force purchased over S 2 billion worth of weapons for 
use in different theaters around the world. The projected need for the quantity of differ- 
ent weapon types is based on an annual Nonnuclear- Weapon Consumables Analysis 
(NCAA) performed by the Directorate of Plans, USAF [Ref. 1]. Unlike other services, 
the USAF relies widely on mathematical programming models in order to optimize the 
allocation of weapons. 

In 1974 RAND developed a nonlinear programming model that optimizes the 
number of different sortie types assigned to several target types by maximizing the mili- 
tary worth of killed targets [Ref. 2: p. 5]. Since each target type was given a different 
target value, the model attempts to assign sorties to maximum value targets first. To 
avoid an undesired concentration of sortie allocations to a few or even one target type, 
a nonlinear objective function was introduced. Within the model only the number of 
available targets and sorties are constrained. The expenditure of weapons is not consid- 
ered. The number of targets one sortie is able to destroy is expressed by an effectiveness 
parameter that depends only on sortie and target type. 

The required input data structure for the RAND-model is a simplification of the 
much more complex data base contained in the Joint Munitions Effectiveness Manual 
(JMEM) used by USAF. The JMEM data base determines effectiveness as a function 
of weather and mission profile (tactic) as well as type of aircraft and type of target. In 
the current operation a model called SELECTOR sorts the JMEM data base so that for 
each sortie-target type combination, all feasible tactics are ordered from the most to the 
least cost-effective, including the cost of aircraft attrition. This list is referred to as the 
Preferred Weapon List. 

The data in the Preferred Weapon List must be reduced to input parameters de- 
pending only on sortie and target type as mentioned earlier. This is basically done by 
selecting the most cost-effective tactic from the list feasible for weather situations con- 
sidered in the model. After the optimization has determined the optimal number of 
sorties assigned to different targets, the number of remaining targets and the expenditure 
of weapons is evaluated. This process is repeated in subsequent time periods with a new 
inventory of sorties and also by recording the remaining number of active targets and 
weapons available. In this way, tactical changes in a given scenario over time are 
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considered by optimizing sequentially for discrete time periods. This process is 
accomplished in one programming model and is called the HEAVY ATTACK model. 
The USAF interest is mainly in the consumption of weapons utilized over all time 
periods. 

The objectives of this Thesis are to include a weapons constraint in the 
RAND-model and to investigate alternatives to the currently used objective function. 
In addition, the RAND-model is expanded so that more available information is in- 
cluded in the optimization in order to gain a higher total military worth of killed targets 
than is currently achieved. Therefore, the consumption of weapons used by less cost- 
effective tactics is investigated when other weapons, used by the most cost-effective 
tactic, are exhausted. As a final consideration, one global optimization over all time pe- 
riods is compared to the current sequential optimization method. Global optimization 
achieves a hisher overall worth of killed targets. However, aainina a higher militarv 
worth of killed targets serves only as an aid in analyzing the predicted need of weapons. 
The value of the revisions suggested in this Thesis have to be measured on their ability 
to satisfy the demands of the USAF and simultaneously meet budget constraints. 
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II. BASIC STRUCTURE OF HEAVY ATTACK 



A. THE ORIGINAL RAND - MODEL 

In 1974 RAND developed a nonlinear programming model whose objective was to 
determine the optimal number of sorties of type i assigned to targets of type j by maxi- 
mizing the total military value of destroyed targets. The relationship between an as- 
signed sortie and a target kill is established by introducing "sortie effectiveness" E t j . The 
parameter E UJ defines the average number of kills that one sortie of type i will achieve 
when it is assigned to targets of type j. 

Definition of index 

i sortie type 

j target type 

Parameter 

Tj total number of type j targets available at the beginning of a time period 

Vj military worth of type j target 

S, total number of type i sorties available 

£N average number of type j targets killed by one type i sortie 



Variables 

SX,j number of type i sorties assigned to type j targets 
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Model 



Max 



-2>, 




SX,J 



S.t. 



X sx u £ 



V i 



fi 



Y, e ‘-j x sx '-i ) - T > 



Vj 



y sx'ij < c x y 5A;j 

j € •/ j 



where J is a subset of all targets of type j and 0 < c < 1. 



0 < SX L 



J 



V i, j 



f,(LE,j*S A’J is a concave function that approaches 1 for large arguments. The 
RAND - model (and HEAVY ATTACK) utilizes a specific analytic from that will be 
examined in detail later. The recipe constraints Y SAV < c x YSAY limit the number 

jsj { 

of sorties of type i which are assigned to a list of targets by a fraction of the total number 
of sorties of type i. Since these constraints are not used by the USAF in their current 
weapon analysis, this inequality will omitted from now on in the Thesis. 
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B. THE ROLE OF SELECTOR 



Based on the information contained in the J.MEM the effectiveness of a sortie de- 
pends on sortie type, target type, weapon type , weather and tactics or mission profile. 

Definition of index 

i sortie type 

j target type 

k weapon type 

w weatherband index 

r index for used tactic 

Definition of parameter 

number of type j targets killed by one type i sortie using tactic r 
in weatherband w 

number of weapons carried by one type i sortie which is assigned to type j 
target in weatherband w and using tactic r 
K.j.r.u, type of weapon which is loaded on sortie i and will be deployed to target j 
by using tactic r in weatherband w 

The JMEM data have too many subscripts to match the required input data structure 
of the RAND - model. The number of subscripts of a sortie needs to be reduced so that 
£, ^depends only on sortie and target type. The first part of the task of reducing the 
number of subscripts from 4 to 2 is accomplished by the sorting program SELECTOR. 
The output data of SELECTOR - referred to as Preferred Weapon List - contains for 
each different sortie - target type combination five distinct items: 

1. The worst weatherband in which a tactic can be used. 

2. The types of weapons that can be allocated. 

3. The relative cost-efficiency of a tactic given by its order on the list. 

4. The number of targets which can be killed by one sortie. 

5. The number of weapons that can be carried by one sortie for each weapon type 
(mixes of weapons are not considered). 
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The data structure of the Preferred Weapon List, which will be used later for the aggre- 
gation of the input data E UJ for the RAND - model, is illustrated by the following ex- 
ample: 

Subset of data from Preferred Weapon List 



i 


j 


r 


w 


IS 


E 




1 


29 


1 


4 


3 


0.137 


4 


1 


29 


2 


3 


1 


0.664 


6 


1 


29 


3 


2 


17 


1.580 


2 


1 


29 


4 


5 


17 


1.600 


2 



Tor example, the most cost-enicient and feasible tactic for weatherband w=3 is 
tactic r=2. Tactic r= 1 is more cost-efficient because it is first on the list, but is only 
feasible in weatherband w = 4 or higher. Weatherband w= 1 expresses best weather while 
weatherband w=6 represents the worst weather. Tactic r= 3 is feasible (a tactic feasible 
in w is always feasible in better weatherbands) but less cost-efficient than tactic r= 2. 

The given data can be represented in the following way: 

Table 1. E U R W - VALUES: Number of targets of type j killed by one sortie of type 



i using tactic r in weatherband w. 



i 


j 


r 


\v= 1 


w = 2 


\v= 3 


\v= 4 


\v= 5 


w = 6 


i 


29 


1 


0 


0 


0 


0.137 


0.137 


0.137 


i 


29 


2 


0 


0 


0.664 


0.664 


0.664 


0.664 


i 


29 


■*> 


0 


1.5S0 


1.580 


1.5S0 


1.580 


1.580 


i 


29 


4 


0 


0 


0 


0 


1 .600 


1 .600 
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Table 2. B, ]fLW - \'ALUES: Number of weapons that are loaded on one sortie of 

type i which is assigned to target type j and using tactic r in weatherband 
w. 



i 


j 


r 


w= 1 


w= 2 


\v= 3 


\v = 4 


\v= 5 


\v= 6 


i 


29 


1 


0 


0 


0 


4 


4 


4 


i 


29 


2 


0 


0 


6 


6 


6 


6 


i 


29 


3 


0 


2 


2 


2 


j- 




i 


29 


4 


0 


0 


0 


0 


2 


2 



Table 3. K ukw - VALUES: Type of weapon that is allocated to a sortie of type i 

which is assigned to a target of type j and using tactic r in weatherband w 



i 


j 


r 


w= 1 


w= 2 


\v = 3 


\v= 4 


w = 5 


\v= 6 


i 


29 


i 


0 


0 


0 


j 


-> 

j 


j 


i 


29 


2 


0 


0 


i 


i 


i 


i 


i 


29 


3 


0 


17 


17 


17 


17 


17 


i 


29 


4 


0 


0 


0 


0 


17 


17 



Since HEAVY ATTACK only considers the tactic at the top of the list for each 
weatherband, and since weapon type is implied by tactics, SELECTOR essentially re- 
duces the number of subscripts from 4 to 3. 

C. DETERMINATION OF E u IN HEAVY ATTACK 

An important assumption for HEAVY ATTACK in order to understand the logic 
behind the aggregation of E tJ is that the weather is not known at the time when sorties 
are assigned to targets. This leads to the condition that the effectiveness of a sortie and 
the consumption of weapons in a particular weatherband has to be proportional to the 
probability that this weather will occur. 

This probability is represented in HEAVY ATTACK by a given distribution of 6 
distinct weatherbands: 

PR w = probability that weatherband w will occur at a certain time in the future, 
w = 1, 2 6. 
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Throughout this Thesis the following distribution is used: 



Table 4. WEATHER DISTRIBUTION IN HEAVY ATTACK: Probability that 

weatherband w occurs when sorties are allocated to targets. 





w= 1 


w= 2 


\y— 3 


w = 4 


\v= 5 


\v= 6 


PR„ 


0 


0.02 


0.14 


0.07 


0.07 


0.70 



Since weatherband w= 1 will never occur, the effectiveness for any sortie in this 
weatherband is irrelevant. It is assumed that any weapon which is feasible for a certain 
sortie - target combination can be used in the weatherband determined by SELECTOR 
or in any better weather (higher weatherband). 

HEAVY ATTACK uses for each weatherband only the top weapon on Preferred 
Weapon List. This means that the model will allocate the most cost-efficient weapon 
feasible in each weatherband. Therefore the data set E iJru can be reduced by the sub- 
script r such that: 

Ejj yW = the effectiveness of the most cost-efficient tactic in weatherband w. 



Table 5. EFFECTIVENESS OF THE MOST COST - EFFICIENT TACTIC: In 

each weatherband w the first effectiveness value in Table 1 greater than 
zero is selected. 





w= 1 


w = 2 


\v= 3 


\v= 4 


\v= 5 


\v= 6 


e;„. 


0 


1.5 SO 


0.664 


0.137 


0.137 


0.137 



Applying the same reasoning on the data set and K, jrw yields : 

D i j w = number of weapons used by the most cost-efficient tactic in weatherband w, 

* 

Kij w = type of weapon used by the most cost-efficient tactic in weatherband w. 
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Table 6. WEAPON LOAD OF THE MOST COST - EFFICIENT TACTIC: In 

each weatherband the first weapon load value in Table 2 greater than zero 
is selected. 





II 

£ 


w= 2 


3 

II 


II 

£ 


\v = 5 


'O 

II 

£ 


b;,,„ 


0 


2 


6 


4 


4 


4 



Table 7. WEAPON TYPE OF THE MOST COST - EFFICIENT TACTIC: In 

each weatherband w the first weapon type in Table 3 not equal to zero is 
selected. 





\Y= 1 


\v= 2 


CO 

II 

£ 


\v= 4 


\Y — 5 


\Y“ 6 


a':,. 


0 


17 


i 


3 


3 





Since each weatherband will occur with the probability PR W . the averaged 
elfectivness must be 



E;j = Yj PR " x £ 'V> = °- 240 
w 

In general the process of obtaining £ (i/ is a little more complicated than described 
above because HEAVY ATTACK is permitted to use tactics lower than first order when 
first order weapon types have been exhausted. This can happen because HEAVY AT- 
TACK is actually a model of protracted war. First order tactics are preferred because 
they represent the most cost-effective tactic. The war may last for several periods (4 in 
this Thesis), and it is possible that certain tactics may not be feasible in later periods on 
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account of weapon exhaustion. Suppose for example, that weapon type 3 has been ex- 
hausted in a previous time period and is therefore no longer available. The top weapon 
for weatherband w=4. 5 or 6 is now weapon type 1. The new effectiveness values 
are: 



Table S. E Uftw - VALUES AFTER WEAPON K=3 IS EXHAUSTED: Number 

of targets of type j killed by one sortie of type i using tactic r in 
weatherband w that is applicable. 



i 


j 


r 


w= 1 


w= 2 


\v= 3 


\v = 4 


\v= 5 


w = 6 


i 


29 


i 


N A 


N A 


N A 


N A 


N A 


N A 


i 


29 


2 


0 


0 


0.664 


0.664 


0.664 


0.664 


i 


29 


3 


0 


1.5S0 


1.5S0 


1.5S0 


1.5S0 


1.5S0 


i 


29 


4 


0 


0 


0 


0 


1.600 


1.600 



Using the most cost-efficient tactic in each weatherband w gives the following ef- 
fectiveness values : 



Table 9. EFFECTIVENESS OF THE NEXT FEASIBLE COST - EFFICIENT 
TACTIC: In each weatherband w the first applicable effectiveness value 

in Table S ereater than zero is selected. 





w= 1 


\v= 2 


w — 3 


\v = 4 


\v= 5 


\v= 6 




0.000 


1.5 SO 


0.664 


0.664 


0.664 


0.664 



which results in the averaged effectiveness: 

E,J = Yj PR " x E U' w = 0,682 • 

VC 

Note that the effectiveness has increased on account of the lack of weapon type 
k = 3! The SELECTOR output is ordered according to cost-effectivness (not effective- 
ness), so it is quite possible that tactics far down in the Preferred Weapon List may ac- 
tually be quite effective. These tactics typically have high associated attrition, but 
attrition is not considered in HEAVY ATTACK once SELECTOR has done its job. 

By considering the same logic, it can be observed that the fourth order tactic on the 
Preferred Weapon List with = 1.600 will never be used. This is because the third 



10 



order tactic uses the same weapon (in this case weapon type k= 17) in at least the same 
worst weatherband as tactic r=4. 

D. TIME IN HEAVY ATTACK 

Once the effectiveness values E t j are evaluated, the required input data is available 
in order to optimize the number of sorties assigned to the different target types. For 
most cases all targets are not killed when the optimization is finished because of the 
constrained number of sorties in the RAND - model. As in a real war scenario, the 
outcome of a given attack will infiuence subsequent target consideration and planning. 
Only the targets that survived the previous attack will be reconsidered. Weapons are not 
resupplied and therefore may become exhausted. The current version of HEAVY AT- 
TACK may actually allocate more weapons in a given period than are available at the 
beginning of the period. This is because there is no explicit constraint on weapon usage. 
The deletion is currently done after each period by computing weapon usage after the 
optimization for the period is finished. However, a weapon will be deleted in the next 
period if it is exhausted at the end of the current period. 

There is no resupply of targets between periods in HEAVY ATTACK, although 
there is a facility for reconstituting targets that have already been killed. This will be 
discussed later. Aircraft are also not resupplied or even directly represented in HEAVY 
ATTACK: the number of sorties available during each period is a direct input. Each 
time period represents an attack which changes the input for the following time period. 

The fact that the importance of a target will change with time is represented in 
HEAVY ATTACK by the option of changing the military worth for each target type 
at the beginning of a new time period. Even though the military worth of a target is 
known in all future periods, the current sequential time optimization only "secs" the 
worth of a target for the current time period. Following from this "myopic" way of 
maximizing the military worth of killed targets it may happen that sorties are assigned 
in a time period to a target type when its military worth is relatively low. A "global" (or 
overall) time optimization can be expected to achieve a higher military worth of killed 
targets. This is discussed later. 



11 



E. THE NONLINEAR MODEL IN HEAVY ATTACK 

The basic structure of tire current model in HEAVY ATTACK for one time period 
is given by: 

Parameter 



Tj number of type j targets available at the beginning of a time period 

Dj number of dead type j targets at the beginning of a time period 

Vj military worth of type j target during the current time period 

Cj target - parameter for type j target 

S, number of type i sorties available for the current time period 

PROP, proportion of S, that can be assigned 

Variables 



SXjj number of type i sorties that are assigned to type j targets 

KILL , number of type j targets killed in the current time period 
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Model 



Max z= 7l}x KILL 



s.t. 



KlLLj = f^Tj, Cj, Dj, x SX U 



V j 



where: 



f( Tj, Cj. Dj, ^ x SA[j 



T - °') * 



(I - 



The above function is the same function as used by RAND [Ref. 2]. 



Vsa; v < PROP, x S, 



V i 



0 < KlLLj < Tj - Dj 



0 < SX U 



vy 
v '• y 



The nonlinear function f(T r c ; , £>„ ]T£ v x .SA' V ) is of the same form as in the RAND 

i 

- model. The number of targets of type j that are killed and the number of sorties of type 
i are constrained. The consumption of weapons is not considered in the model itself. 
After the optimal numbers of sorties are determined by the optimization, the consump- 
tion of the different weapon types is evaluated by: 
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{ consumption of weapon}^ = / t / SXq x ( £^PR W x B lt j >w j 

i j s w / 

where the sum is over all { i, j. w } such that k = K' J W . 

F. TARGET RECONSTITUTION IN HEAVY ATTACK 

The ability to reconstitute killed targets is a common fact in a modern war. HEAVY 
ATTACK records the number and type of targets as well as the time period when they 
are destroyed. After each optimization, it determines if targets can be reconstituted and 
evaluates the maximal number that are possible. A major task in this Thesis has been 
to determine the conditions under which reconstitution is allowed to happen by analyz- 
ing the responsible part of the HEAVY ATTACK source code. HEAVY ATTACK'S 
logic seems to be as outlined below: 

Definition of index 



j 


target type index 


V j 


p- pp 


time period index V p. pp e (1. 2,. 


•• , «} 


Parameter 






TIME p 


length of time period p in days 


V p 


RECOX j 


minimum number of days a target has to stay dead 


v; 


QTYj 


maximum number of targets j that can be reconstituted in 30 days 


vi 


Aggregated parameter 




PERUPjp 


index of the last time period considered for reconstitution. 
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If a target of type j is killed in time period PERUP or earlier, then there is sufficient 
time available to reconstitute the target so that it once again will be available in period 
p+1. The parameters TIME, and RECO.X, determine PERUP, according to the fol- 
lowing formula in HEAVY ATTACK: 

p + 1 

if RECOXy < Yj TlME p . - CEIL (0.5 x TIME-) V j, p < p < n 
otherwise p '=p 

where the function CEIL rounds a real number to the next higher integer value. 

L ; - iP indicates whether targets killed in period p are eligible for reconstitution in 
period p and therefore: 

p 

PERUPjp = Y k j,p, P v J- P < n 

p=i 

Note that always PERUPj p < p. 

Variables 



Let 




KILLjp number of targets type j killed in time period p V j, p 

REBUILDj p maximum number of targets of type j that are reconstituted 

as live targets in time period p+1 V j, p < n 

Conditions for Reconstitution 

A killed target of type j can be reconstituted if the following 4 conditions are true: 

1. at least a fraction of target j was destroyed in the previous or the current time pe- 
riod p, 

2. it has been dead for more than some defined time, 
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3. the total number of targets being reconstituted has to be less than the total number 
of targets which exceeds the minimum dead time 

p per L'P ij, 

Y REBUILD jp , < y KILLjp, V j, p < n 

p'=\ p'=\ 



4. the maximum number of targets type j which can be reconstituted at the end of 
each time period p is given by: 



QTYj 

REBUILD j p < -pp x TIME p+] 



where 



OTYj 

30 



represents the reconstitution rate per day. 



V j, p < n 



This leads to the following submodel: 



max 



Z Yj rebu,ld u 

j p 



s.t. 



p PERUP„ 

y^REBUILDj p , < y KILLjp. V j, p < n (A) 

p '= i p '= i 

QTYj 

REBUILD jp < x TIME p+] V p < n (B) 

The interpretation of (A) is that the number of targets of type j rebuilt in period p or 
before cannot exceed the total number of targets that are killed during or before period 
PERl'Pjj . The interpretation of(B) is that the number of targets of type j rebuilt in 
period p cannot exceed a certain quantity depending on the length of period p and on 
the target type. There are no targets reconstituted in the last time period p= n . 
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III. BOUNDS ON WEAPON CONSUMPTION 



A. INTRODUCTION OF A WEAPON CONSTRAINT 

A desired improvement for the current HEAVY ATTACK model is to add an addi- 
tional constraint on the utilization of weapons inside the RAND - model. 

Two important facts should be recalled: 

1. For each sortie - target combination { i, j } and each weatherband there is at most 
one weapon which can be used. 

2. Averaging over all weatherbands is related to the probability that weatherband w 
might occur at the time sortie type i is assigned to target type j. 

Let the upper bound on weapon consumption be defined as: 

ITPfc total number of weapons of type k available 
The required constraint for the consumption on weapons is then: 

E E SA 'u * ( E™- x K].~ ) z n -r k v k 

i J \ * / 

where the sum is over all { i, j. w } such that k = A'Y * 



B. REVISED MODEL OF HEAVY ATTACK 

Reconstitution can be included in the RAND - model. Instead of considering re- 
constitution as a computational "bookkeeping" process, it can be part of the optimiza- 
tion. To accomplish this, it is necessary to define a new variable for the number of dead 
targets such that the time period as an additional dimension is represented by a second 
subscript: 

D jr is the total number of targets of type j killed in time periods < p less the number 
of targets that are reconstituted during this time : 



p - 1 

Djj, = V (KILLjp- - REBUILD »• p .) 

p'= i 



V j. P 
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The military worth of a target is also time dependent: 

1', p military worth of a target type j in time period p 

Embellished Thesis Model (solved sequentially for p = 1,2, 3,..., n) 

Max z p = Tj V j.P x KILL J,p) 

J 

s.t. 

KILLjp = /|t ; , Cj , D jp , ^SXjj x {Yj PR w x Kj, 
where : /(...} is one of three functions discussed in the next chapter. 




V j 



KILLjp < Tj - Djp 



p - 1 



Dj <p = y (KILL ] p , - REBUILD^) 






PER UP , 



y^REBUILIhp. < V KILLjp, 
p '= i />'=i 



YsA’-j < TAO/ 5 ,- x S,- 



p K x Aj.J f ^ wp k 



zeKj * & 

i j (. V w 

where the sum is over all { /, j , w } such that k = A'*,- w 



V j 



Vi 



Vi 



V i 



V A 



IS 



0 < SA'ij 




v /, J 


0 < KlLLjp 




V j 


0 £ Dj, 




vy 


0 £ REBUILD, p 




* j 


where the upper bound on 


REBUILD] p is such that: 




r< QTr, 

REBUILD] p 30 

1=0 


X 

+ 

A 

D 


vy 


if p = n 


v j 



The model was written in the General Algebraic Modeling System f GAMSj [Ref. 3]. All 
optimization problems throughout the Thesis are solved with the nonlinear programing 
solver Ml SOS - Version 5.0 [Ref. 4], A database for 2 sortie-. 26 target- and 29 
weapon-types was provided [Ref 5] in order to compare the results by using three dif- 
ferent objective functions, each over four time periods. 
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IV. LINEAR VERSUS NONLINEAR MODEL 



In this chapter the derivation of the nonlinear objective function used by RAND is 
given. In addition two alternatives are represented by introducing the Washburn- 
Equation and the linear case in which the number of killed targets is proportional to the 
number of assigned sorties. Each of the three objective functions is used in the model 
described in the previous chapter for sequential}' optimizing sortie assignments over four 
time periods. In order to compare the effect of the three objective functions, a meas- 
urement for the diversity of the allocated kill capability is defined. 



A. RAND EQUATION 

If K represents the total number of killed targets of type j then the objective function 
used in the RAND - model can be derived from the differential equation: 



4 % 

d Xj 



= 1 — C, X 



T J 



(A) 



where Xj 



'y' * E-q x SXq and 0 < Cj < 1 



The differential equation (A) with the initial condition A'(A' = 0) = D, has the solution: 

T: C D, _± y ) 

K ] = — x | 1 - (1 - cj x — ) x e Tj A 

Instead of bounding A' by 



Dj < Kj < Tj 

let KILLj be the number of targets killed in excess of D, : 

KlLLj = Kj - Dj 



so that 
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which leads to the final result: 



0 < KILLj <Tj - Dj 



KILL, - ( - D-'j x ( 1 - e- T * 

B. LINEAR EQUATION 

A special case for the differential equation (A) appears when c, 
then 



d Kj 
~dXj 



1 



which yields: 



so that 



or by using 



so that 



Ay = Ay + Dj 
Dj < Kj < Tj 

KILLj = Kj - Dj 



0 < KILLj < Tj — Dj 

where the final solution represents the linear case: 

KILLj = Xj 
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Figure 1 illustrates the influence of the target parameter 
Cj on the function KILL - t = /(A'). 



KILL CAPAPILITY = X 

NUMBER OF KILLS GAINED = KILL 




KILL CAPABILITY 



Figure 1. Influence of the target parameter c on the RAND-Equation: The sol- 

ution of the differential equation used in the RAND-model is graphically 
shown for 4 different target parameters c. 



The parameter c J has no direct physical motivation. The model considered in the next 
section also contains a single parameter, but the parameter can be motivated physically. 



C. WASHBURN EQUATION 

d K 

The Washburn - Equation [Ref. 6: p. 25] defines the differential — — in the fol- 
. . d A 

lowine wav: 



d Kj 



= Probability { attacking a live target } 



or equivalently: 



d K] { number of live targets } 

d Xj { number of targets that look alive } 



This leads to the differential equation: 



dKj_ 

d Xj 



Tj - Kj 



Tj - Kj + c/.j x Kj 



(B) 



where a, is a constant proportion of killed targets, which have the property to appear live 
to a potential attacker. 

The differential equation (B) with the initial condition A'(A' = 0) = D, has the sol- 
ution: 



K J = T J * 1 1 



3 . 

Tj 



(1 - 2 ,)x (K, - D.) - Xj 



X e 



y j x T j 



Using KILLj instead of A' such that: 



KILL] = Kj - Dj 



leads to the implicit solution for the Washburn - Equation as: 

/ (1 -a,)x KILLf-X, 

KILL: — (T, - D,) x ( 1 - e TJTt, 



The difference between the two differential equations (A) and (B) for two different target 
parameters is shown in Figure 2 on page 24 . Observe that for target parameter c close 
to 0 or 1 the Washburn-equation tends to behave similarly to the RAND-equation. 
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Target parameter a is denoted in the figure by c. 




Figure 2. RAND- and Washburn-Diff. Equation with varied parameter c: The two 

differential equations are shown for 2 different target parameters c. Be- 
cause the solution of the Washburn-Equation can be given only in an 
implicit form, the differential equations are shown rather than their sol- 
utions. 



The influence of the three different objective functions on the RAND-model using 
the same input data is shown in Figure 3. 

The total worth of killed targets decreases with time for each objective function. 
The main reason for this is that in the first time period sorties are assigned to those 
target types for which the effectiveness is highest. When all targets are killed, sorties are 
then assigned in the following time periods to the remaining targets for which the effec- 
tiveness is less. As a result, more and more sorties need to be allocated in order to gain 
the same number of killed targets. The number of reconstituted targets available at the 
beginning of the second or third period is relatively small or even zero and can therefore 
be neglected at this point. Since the variation in the number of sorties and in the mag- 
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nitude of the target values is too small to compensate for this effect, a declining trend 
in the objective function value over time for all three cases is observed. 

Note that the Washburn-Equation always yields a smaller value than the 
RAND-Equation. This follows from the fact that the Washburn-Equation declines faster 
than the RAND-Equation for the same target parameter c as shown in Figure 2. The 
linear equation is larger than either one. The most important difference is not in the 
absolute level of target value killed, but rather in the influence of the objective function 
on the distribution of sorties over targets. This subject is taken up in the next section. 



Objective Function Type: 

R = RAND 
L = Linear 

W = Washburn 




F'g ure 3. Total Military Worth of Killed Targets: represented for each different 

objective function and each time period by the height of the respective 
block in the figure. 
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D. DIVERSITY OF KILLED TARGETS 

An important reason for USAF to use a nonlinear objective function is to avoid an 
undesired concentration of attacking sorties on a few targets. In analysing the effect of 
the three different objective functions on the optimization, a measurement is needed in 
order to indicate how many of the allocated sorties are spread over different targets. 

In information theory the function 

%) = (pt x log j- ) 



where p = (/>,, p 2 , ... , p n ) and ^p, = 1 

i 

is used to express the diversity or "entropy" of the probability distribution p = {p, }. 
Observe that /i(p) = 0 when p concentrates all probability in one element. The maxi- 
mum possible value when p has n elements occurs when they are all equal, in which case 
/j(p) = log n . The diversity h(\) of an arbitary set { x ; } of nonnegative members can 
be measured by simply normalizing them so that they sum to 1 and then computing 
entropy: 



y Xj x log 

j 

h{\) = 

j 

The diversity of values h{\) gained from the same input data and model as used in the 
previous chapter is depicted in Figure 4. Since the number of targets n equals 26, the 
maximum diversity value will be 





/j (*) ma x = 3.26 

Figure 4 makes it clear that the Linear objective function has a lower diversity value 
than the other two. This is to be expected, and in fact one of the main reasons for using 
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a nonlinear objective in the first place was to avoid low diversity values. However, note 
that: 

1. The Linear diversity is not 0; that is, several target types are still attacked. 

2. None of the objective functions achieves complete (3.26) diversity. 

The differences emerge most strongly in period 3. Only 4 target types are attacked when 
the linear model is used, or 6 with the RAND-model. 16 different target types are at- 
tacked when the Washburn-equation is used; this is in keeping with the idea that the 
Washburn-equation is the most "non-linear" of the three (see Figure 2). The three 
models differ much less in period 1.2 or 4. 
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Objective Function Type: 

R = RAND 
L = Linear 

W = Washburn 




Figure 4. Diversity of killed targets for different objective functions: The height 

of each block illustrates to how many different target types (out of 26 ) 
sorties are allocated at different time periods by using each of the three 
objective functions. 
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V. ALLOCATION OF SECONDARY WEAPONS 



A. COST-EFFICIENCY VERSUS KILL-EFFECTIVENESS 

Cost considerations are finished once SELECTOR has established the Preferred 
Weapon List. Although this list contains different tactics, ordered in terms of cost- 
efficiency, HEAVY ATTACK only uses the top one on the list which is feasible. The 
only time at which HEAVY ATTACK may proceed to a succeeding tactic appears, as 
mentioned before, when a weapon has been exhausted in earlier periods. 

As a second revision of HEAVY ATTACK, the model is changed to continue target 
attacks after the weapon type used by the most cost-effective tactic has been exhausted, 
using those weapons still on hand. 



B. A NONCONVEX CONSTRAINT 

The model discussed in the previous chapter requires that only the tactic on the top 
of SELECTOR'S Preferred Weapon List can be used. Once the corresponding weapon 
type is depleted further attacks by that sortie type in that weatherband against that tar- 
get type are impossible. The idea in this section is to relax this strict requirement to 
permit using whatever tactic is highest on SELECTOR'S list among those whose weapons 
have not been exhausted. 

Implementing this logic in the existing model requires a modification of the variable 



SX,^. r w = number of sorties of type i assigned to target of type j which use 
tactic r in weatherband w 

The probability that all sorties of type i assigned to target of type j will attack the target 
in weatherband w has to be equal to the probability that weatherband w occurs at that 
time: 



S V 
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L'pon these redefined variables for the number of assigned sorties, it is possible to de- 
termine the utilization of each weapon type: 

let ]VEAP k be the consumption of all weapons of type k 



where the sum is over all { /, j, r, w } such that k = Kij tr w . 

In order to assign sorties using less cost-effective tactics, SX iJ%riM must be 0 unless the 
weapon types corresponding to all more cost-effective tactics are exhausted. The fol- 
lowing constraint will enforce this logic: 



The above constraint requires that at least one of the two factors on the right hand side 
of the equation equals zero, so either no sorties are assigned (first factor zero) or else 
all more cost-effective weapons are exhausted (second factor zero). The constraint thus 
enforces the desired logic, but there is a disadvantage in using it. The disadvantage is 
that the function on the right hand side of (C) is not only nonlinear (products of vari- 
ables are involved) but nonconvex. Without constraint convexity, there is no guarantee 
that the locally optimal solutions achieved by the MI.XOS solver are globally optimal. 
There is some evidence, however, that globally optimal solutions are actually being at- 
tained. For one thing, employing constraint (C) always results in a higher objective 
function value than when only the most cost-efficient tactic is permitted. In addition, 
some experiments were performed where the improved model was changed into a linear 
model by linearizing the objective function at the optimal solution. The nonconvex 
constraint was then converted into a linear constraint by using integer variables. The 
optimal solution of this linearized model was identical to the solution gained by the 
nonlinear model with the nonconvex constraint. 



then 




V k 



r— 1 




where k = K iJr , <w . 
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C. REVISED MODEL 

The mathematical model is solved sequentially for p = 1. 2 ,..., n. 



Max 2 , = x K1LL hp) 



s.t. 



KILL u = 



7~ ~ D Ip ) x (l - e Tj x */) 



V j 



where A} = V Z S x SX UsJ 



KILL hp < Tj — D hp 



V j 



P -\ 



Dj P = V (KILLj p, - REBUILDj p .) 

/>'=! 



Vy 



PER VP, 



'y^REBUILDj p , < y KILLjp, 



P'= 1 



P'=l 



V j 



EEE sa w £ x s < 



V i 



I! EAR k — ZZZZ x SXijs,w) 



‘ j 



V A 



where the sum is over all { /, J, r, tv } such that A = M 
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